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As a first step towards a microscopic understanding of the effective interaction between col- 
loidal particles suspended in a solvent we study the wetting behavior of one-component fluids 
at spheres and fibers. We describe these phenomena within density functional theory which 
keeps track of the microscopic interaction potentials governing these systems. In particular we 
properly take into account the power-law decay of both the fluid-fluid interaction potentials 
and the substrate potentials. The thicknesses of the wetting films as a function of temperature 
and chemical potential as well as the wetting phase diagrams are determined by minimizing 
an effective interface potential which we obtain by applying a sharp-kink approximation to the 
density functional. We compare our results with previous approaches to this problem. 
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I. Introduction 



If colloidal particles are dissolved in a solvent consisting of a binary liquid mixture, 
necessarily one of these two components is preferentially adsorbed on the spherical surfaces 
of the colloidal particles. Near a first-order phase separation of the bulk solvent this 
adsorption may lead to the coating of the colloids by wetting films which snap into a 
bridgelike structure if two colloids come close to each other. This coagulation can become 
so pronounced that it results in flocculation . Such experiments have inspired several 
theoretical investigations [|^-[T6[| which have shed light on various features of these complex 



phenomena. 

Although there are indications that for higher concentrations of colloidal particles their 
effective mutual interactions are beyond a superposition of effective pair potentials |TT[, for 



low concentrations the knowledge of the effective interaction between two isolated colloidal 
particles at a fixed distance immersed into the solvent is important. Moreover, this latter 
configuration is interesting in its own, because the bridge formation between individual 
pairs is not solely accessible indirectly via flocculation but can be studied also directly 
with atomic or surface force apparatuses [T^. Besides colloidal systems these problems 
are in addition relevant for technical applications such as the contact between lubricant 



films on the disk surfaces in magnetic recording with the recording head [0, sintering 
|T9| , latex film formation catalyst wetting efficiency in trickle-bed reactors |^, or 



for the bonding mechanism of emulsified adhesives in granular and fibrous substrates 

As a prerequisite for understanding the effective interaction between such spherical par- 
ticles immersed into a solvent one must know the structural and thermal properties of 
these systems in the case that the distance between two colloidal particles is macroscopi- 
cally large. Near a first-order phase transition of the solvent this amounts to studying the 
wetting behavior of the solvent at the surface of a spherical substrate. The paradigmatic 
case is the formation of a liquidlike film at the curved substrate-vapor interface upon ap- 
proaching the liquid- vapor coexistence curve in the bulk of a one-component fluid from the 
vapor side. This problem has been studied theoretically by several authors [^3| -p9| based 



on various versions of phenomenological models. They confirm the general expectation 
(see Subsect. X.B in Ref. [Q) that the positive curvature of the surface of a substrate 



prevents the build-up of macroscopically thick wetting films because, in contrast to a pla- 
nar geometry, the area of the emerging liquid-vapor interface increases with the thickness 
of the wetting film. This increasing cost of the free energy of the film has effectively 
the same consequence as if the bulk fluid is kept off liquid-vapor coexistence. Therefore 
continuous wetting transitions are eliminated and first-order wetting transitions are re- 
duced to quasi-first-order transitions between small and large, but finite, film thicknesses, 
smeared out due to the finite size of the substrate area. 

However, beyond these general aspects it is important to know to which extent the film 
thicknesses of these wetting films are limited and how these limitations depend on the 
radius of the colloidal particles, on the character and the form of the substrate potential 
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and the interaction potential between the solvent molecules, and on the size of the solvent 
molecules. The aforementioned phenomenological models do not allow one to answer 
these questions because they do not keep track of the microscopic details of the system. 
In particular the experience with wetting phenomena on flat substrates tells that it is 
essential to take into account equally the power-law decay of the substrate potential and 
of the interaction potential between the fluid particles pO| , |3T| . Inter alia, the long range 
of the forces between the fluid particles causes the breakdown of a gradient expansion for 
treating the deviations of the emerging liquid-vapor interface from its flat configuration 
and leads to a nonlocal Hamiltonian |^2|]. For simple geometric shapes of the interfaces 
involved this long-ranged character of the dispersion forces acting in and on fluids is taken 
into account by the Dzyaloshinskii-Lifshitz-Pitaevskii (DLP) theory |^ which has been 
applied to spheres and cylinders |Q. For a given system this allows one to compute 
the cost in free energy Qs{l) to maintain a thick liquid film of prescribed thickness / 
adsorbed on the substrate in terms of the frequency dependences of the permittivities 
of the substrate, the bulk liquid, and the bulk vapor. Although this approach has the 
advantage to take into account many-body forces and retardation, it suffers also from some 
shortcomings, (i) The DLP theory gives access only to the leading asymptotic behavior 
of Qs{l —>■ oo) and thus does not allow one to describe critical wetting transitions which 
result from the competition between the leading and next-to- leading order term in Qs{l) 
for large / [^,0. (ii) Details of the substrate potential, which are important for first- 



and second-order wetting transitions pO| , |3T| , are not captured by the DLP theory, (iii) 



exhibits an unphysical divergence for I because this approach ignores the 
repulsive part of both the substrate potential and the fluid-fluid interaction potential and 
thus the structure of the emerging substrate-liquid interface, (iv) The broadening of the 
emerging liquid-vapor interface upon raising the temperature cannot be accounted for. 
(v) The dependence of Qi^^^\l) on temperature and chemical pressure, in particular 
close to phase transitions in the bulk fluid, is not transparent, (vi) In the present context, 
as expounded above, the most severe drawback of the DLP approach is, that the DLP 
results for wetting of a single sphere or cyhnder cannot be used as a building block 
for investigating later on the effective interactions between two such objects because the 
shape of the bridging wetting film between them is not known in advance which precludes 
the practical application of the DLP theory to the problem described in the beginning. 

Therefore in this paper we compute (Sect. II) and discuss (Sect. Ill) the effective in- 
terface potential f2s(/) for the wetting of a single sphere and of a single cylinder on the 
basis of density functional theory [^. This microscopic approach allows one to address 
the aforementioned points (i)-(vi) and overcomes the shortcomings of the phenomenolog- 
ical theories. The price to be paid is that the results do not account for the effects of 
many-body forces. 

In the beginning of this introduction we have focused on the physical interest in spheres 
exposed to a fluid. It turns out that there is also substantial interest in the cylindrical 
geometry. In view of its importance for technical processes such as the lubrication of 



textile fibers, optical fiber processing, and the formation of fiber-reenforced resins [BB 
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there are many studies of the statics and dynamics of adsorption on cyhnders [p|, p6| , p7| -|47[ . 
In addition the adsorption on cyhnders is important for the vibrating- wire microbalance 
as a standard technique to study wetting phenomena However, one has to keep in 
mind that the convoluted surface of graphite fibers impedes the interpretation of these 
experiments on the basis of theoretical models which assume a smooth cylindrical surface 
(Such surface inhomogeneities may even play a role for spherical colloidal particles 
For the same reasons given above for the spherical geometry we apply a microscopic 
density functional theory also to the problem of wetting of wires in order to obtain a 
detailed picture of their behavior in systems with long-ranged interactions. 

We conclude the introduction with two remarks. First, the present study is dealing with 
volatile liquids, i. e., it is important that the liquid phase and the surrounding vapor are 
in thermal equilibrium. Second, we do not address the problem of interfacial wetting at 
curved surfaces, i. e., the substrate is passive. Therefore our results do not apply directly 
to such phenomena as surface melting of small particles and wires. 



II. Effective interface potential from density 
functional theory 



A. Density functional 



As for the description of wetting phenomena in planar geometries we describe the one- 



component fluid by a simple grand canonical density functional |30,35 



^][{p(r)};T,/i]= / £rU[{pir)},T] 



+ 



- f d^r f d^r'w {\r - r'\) p{r)p{r') + f d^rp{v){p^v{v)-p) 
2 Jv Jv Jv 



(2.1) 



is a macroscopic but flnite volume into which the substrate S is embedded but is 
not a part of it. [{p(r)}; T, /i] is the grand canonical free energy for a given density 
configuration p(r) as a function of the temperature T and the chemical potential p. We 
assume a spherically symmetric pair interaction potential w{r) for the fluid particles 
which according to standard procedures [^] is divided into a short-ranged repulsive and 
a long-ranged attractive part w(r) with w{r oo) ~ r~^. Here we do not consider 
retardation effects which would lead to w{r oo) ~ r^''. The divergent repulsive part 
of the interaction potential gives rise to a reference free energy which is mapped onto 
fh{p,T) as the free energy density of a homogeneous system of hard spheres with an 
effective diameter d{T) |^ and number density p. An appropriate expression for fh is 
given by the Carnahan-Starling formula 
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Mp, T) = ksTp (iniv) - 1 + ^2.2) 

where t] = ^pd{TY is the dimensionless packing fraction. The second term in Eq. (2.1) 
accounts for the long-ranged attractive interaction which we treat according to the Barker- 



Henderson scheme [51 



w{r-) = 9(r — a)w{r), (2.3) 

where 9 denotes the Heaviside function. For many purposes it may be useful to adopt 
the Lennard- Jones potential 0Lj(r) as the actual full interaction potential w{r): 

Mr)=0Ur)=46((^)''-(^)') (2.4) 

with depth e and diameter a of the fluid particles so that d{T) = dr (l — e"'^''^^-^^'")). 
For this interaction potential Eq. (2.1) leads to a phase diagram of the bulk fluid with 



a critical point at keTje = 1.09858, p^a^ = 0.27201, and pje = -3.77713 H. The 
fluid is exposed to a substrate which in a first approximation one may consider to follow 



from a superposition of Lennard- Jones interaction potentials 0?; between wall and fluid 



particles. 



^^i-^e./frn'-f^) ). (2.5) 



i. e., v{r) = d^r'(f)^j{\r—r'\) where S is the substrate volume. Ignoring the exponentially 
decaying corrugation effects of an atomically structured substrate, f (r) depends only on 
the radial distance to the substrate surface, v{r) = v{r). In general, however, v{r) is 
more complicated than a linear superposition of Eq. (2.5) and enters the problem as a 
free parameter function. The equilibrium number density, which minimizes Eq. (2.1), 
depends also only on r and yields the corresponding grand canonical potential of the 
system. 

For the above interaction potentials this minimum can only be determined numerically. 
Here, however, we construct an effective interface potential fls{^) which is the cost in 
free energy to maintain a liquidlike film of a prescribed thickness / near the substrate 
although the bulk fiuid is in the vapor phase. The equilibrium thickness /q minimizes 
this effective interface potential. This approach yields a transparent insight into the 
various wetting transitions which can occur ||5^. It has turned out that the so-called 



sharp-kink approximation, in which the trial-functions p(r) compatible with the thickness 
I are approximated by piecewise constant density profiles, yields a surprisingly accurate 
expression for the actual effective interface potential fis(/) (In the present context 

there are also approaches beyond the sharp- kink approximation These numerical 

studies yield valuable information about specific systems but no overall picture for the 
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wetting phenomena near curved surfaces.) For spherical and cyhndrical substrates the 
sharp- kink approximation amounts to using the following trial function (see Fig. 1(a)): 

p(r) = [piQih - r) + pgQir - h)] Q{r - n). (2.6) 

h and tq denote the radius of the emerging gas-liquid interface and of the substrate, 
respectively, so that / = h — tq > dw The repulsive part of the substrate potential 
leads to an excluded volume close to the substrate. This is taken into account by the 
radius ri = ro + dw'^ approximately d^ = |(o" + <^wf)- Actually dy^ is given by the zeroth 
moment of the wall-hquid interface profile Pwi{z), d^, = dz{l — pwi{z)/pi), at gas- 
liquid coexistence p = po{T) [55(c)]. pi = pi(T,p = po(T)) denotes the equilibrium hquid 
density at two-phase coexistence, whereas Pg = Pg(T,p) is the actual vapor density. From 
Eq. (2.6) one obtains the sharp-kink profile for the flat geometry if one sets ro = and 
replaces r by the distance z perpendicular to a half-space filled with substrate atoms (see 
Fig. 1(b)). 

If Eq. (2.6) is inserted into Eq. (2.1) Q can be decomposed into a bulk contribution 
proportional to V^^\ i = s,c, and into subdominant terms. The superscripts s and c 
correspond to a sphere and a cylinder, respectively, so that V^'^^ = ^{L'^ — rj]) and 
V^(c) = 7rM(L^ — Tq); M is the length of the cylinder. In the cylindrical case we consider 
the thermodynamic limit M ^ oo which leaves one with a wire. In the planar geometry, 
i = p, V^P'' = AL where A is the substrate area. In all cases L denotes a macroscopic 
system size (see Fig. 1). The systematic application of the sharp-kink approximation to 
Eq. (2.1) leads for all geometries i = s,c,p and independent of the form of w{r) and v{r) 
to (sometimes we omit the superscripts) 

fi[p(r); T, p] = VQhiPg, T, p) + S,,,,, + 0,(Z). (2.7) 

is the bulk free energy density of the vapor phase. The finite volume V^"^^ introduces 
an artificial gas- vacuum interface which gives rise to the surface contribution Tig^^ac- This 
and other surface contributions Ea^^(r) are proportional to the surface area A[r) of the 
corresponding interface with radius r so that Wmr^oo^a,^^^) / A{r) = aa,i3 is finite and 
is the surface tension between phases (substrate, liquid, or vapor) a and /3 in planar 
geometry: 

S„,^(r) = A(r)a,,^(r). (2.8) 

We denote the r-dependent surface tension as aa,i3{r). The third contribution in Eq. (2.7) 
is the effective interface potential Cls{l) which contains the dependence on the film thick- 
ness /: 

n,{l) = VugAQh + ^Lg + ^Ls + oo{l) (2.9) 

In Eq. (|2l9|) VugAQb = Vuq{pi - Pg){po - p) = VugApAp is the cost in free energy to fill 
the volume V^-^-* = ^{h^ — f^o) and VJ-^'' = 7rM(/i^ — ''^o); respectively, with liquid although 
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/i and T favor the vapor phase. Upon approaching hquid- vapor coexistence fio{T) from 
below, i. e. fi < /iq, AQ), is positive, so that VugAQi, increases with the radius of the 
hquid volume. Thus for all geometries the formation of an macroscopically thick layer is 
prohibited as long as one is off coexistence. 

If the equilibrium value /o(T, /i), which minimizes fls{l), diverges for /i — > fio(T) one 



has a so-called complete wetting transition [^0|,^. However, for curved geometries a 
macroscopically thick film cannot emerge even when coexistence is reached because the 
free energy contribution due to the liquid-gas interface, E/^g = A{h)ai^g{h), increases with 
the area A{h) of the interface, i. e. with the film thickness I. Therefore the equilibrium 
film thickness Iq on spherical and cylindrical substrates will always be limited to a finite 
value [Q. The free energy of the liquid-solid interface does, by contrast, not depend on 
/, so that s = v4(ro)o"i^s(ro) is a constant with respect to the minimization of ^^(0- 

As long as the liquid-gas and the solid-liquid interface are not macroscopically far apart 
from each other, there will be a nonzero interaction between them. This is taken into 
account by the last term in the effective interface potential, lD{1), which by construction 
vanishes for / ^ oo. In the case of a planar substrate, interfacial phase transitions along 
the coexistence curve fio{T) are determined by uj{l) only pO| , |3T| , |54| . For spheres and 
cylinders, however, the wetting behavior at A/i = yUg ~ = is governed by as well 

as Elg. 

It is convenient to normalize the effective interface potential Cls{l) by the area of the 
substrate surface, A{rQ), so that (see Fig. 1) 

^^''^ AQ, + ^^^^ai,g{ro + l)+aUro)+uj{l) (2.10) 



A{ro) " A(ro) 
where 

^(0 = J^^i^)- (2-11) 
The geometrical prefactors Viiq/A{rQ) and A{ro+l)/A{rQ) have the following explicit forms: 

Vii^ _ ^0 / /i + iy^'_i j (2.12) 



A{ro) T+1 \ \ ro 

and 

A{ro + 



:2.13) 



with r = 2 for a sphere and r = 1 for a cylinder. Before determining the surface phase 
diagrams from the minimization of i7s(/) we first discuss the two /-dependent terms g 
and 00 (l), which are interesting in their own right and allow for a comparison with results 
obtained by other approaches. 
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B. Liquid-vapor and solid-liquid surface tension 



From g we obtain an expression for the curvature dependence of the hquid-vapor 
surface tension of a spherical drop and of a liquid thread. Whereas apart from Ref. [p7 



we are not aware of studies of liquid fibers, the surface tension of liquid drops has been 
subject to partially controversial discussions ^E\. Here we calculate the surface tension 



within our approach for both cases explicitly in order to work out the influence of the 
curvature. 

Within the sharp-kink approximation of Eq. (2.1) the planar liquid- vapor surface ten- 
sion is given by |0 



sS = -7T(Apr/ dzt(^\z), (2.14) 



'>9 A ''^ 2 

where t^"^ (z) is the interaction potential of a fluid particle at distance 2; > of a half-space 
filled with fluid particles, too, interacting with a pair potential w: 

&\z) = I dz' I d\ w (y/rj + ^'2^ , (2.15) 



'R2 

= x^+y^. If one inserts the intermolecular potential w chosen above, one obtains for the 
surface tension of the planar interface aj^^ = |7re(T^(Ap)2. Close to this expression for 
ai^g vanishes ~ t where t = {Tc — T)/Tc — >■ instead of the expected mean- field behavior 
because the sharp-kink approximation fails to take into account the broadening of the 
interface on the scale of the bulk correlation length ^ ~ t^^^^. Thus Eq. (2.14) takes into 
account the long-ranged character of the forces but should not be used close to Tc. The 
corresponding expressions for a cylinder and a sphere of radius h have a similar form: 

and 

t^^\r; h) is the interaction potential at distance r > h from the axis of an infinite liquid 
cylinder of radius h: 

/+00 i'2w i-h 

dz j d(j) r' dr'w ( ^r^ + r'^ - 2rr' cos(j) + zA . (2.18) 
00 Jo Jo ^ ' 

Correspondingly, t^^^{r\ h) is the interaction potential at distance r > h from the center 
of a liquid drop of radius h: 
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o ph pr+r' 

t('){r;h) = — dr' r' ds sw{s). (2.19) 

r Jq Jr~r' 

If one inserts Eqs. (2.3) and (2.4) into Eq. (2.19), one obtains for r ^ h 



3 (r/a)^ \r 

i. e., viewed from a large distance the liquid drop appears as a single molecule with 
an effective radius h, exerting a van-der-Waals interaction on any other fluid particle. 
Equations (2.16) and (2.17) are valid for a general form of the interaction potential w{r). 
For the specific choice given by Eqs. (2.3) and (2.4) the expansion in powers of the 
curvature 1/h yields the following expressions for the surface tensions (see Eqs. (A7) and 
(B35) in Appendices A and B, respectively): 

= (2.21) 
^S('» - (^-l'm-^C>2,±)-0(^)), ,2.22) 



6 (h/a)^ {h/af V3 144; \{hla 



and 



If one compares these results with the usual ansatz for the surface tension of a drop 



-SW = -g(l-^^ + S + ---)^ (2-24) 



25 a 

first, one notices that our approach predicts that the surface tensions of the curved inter- 
faces are smaller than the planar one and that both attain the latter one asymptotically, 
\m\h^oo (^[^f^ih) = ai^g . Second, since our density profile (Eq.(2.6)) is antisymmetric 
around r = h and does not contain a symmetric component, within our approach the 



Tolman length S must be zero p8[. Third, the long-ranged forces give rise to logarithmic 



corrections which dominate the usual subdominant term ~ ah ^. Within our approach 
al^g is invariant with respect to interchanging liquid with vapor and therefore they do 

not depend on the sign of curvature of the interfaces. The actual behavior of <y\^g{h) and 

<j\^g{h) as obtained from Eqs. (2.16) and (2.17) is shown in Fig. 2. 

The solid-liquid surface tension entering Eq. (2.10) depends only on tq and thus it is 
constant with respect to the minimization of VLs{l). Within the sharp-kink approximation 
(Eq. (2.6)) it has the following form (r = 2 for a sphere, r = 1 for a cylinder): 
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-^MV^::l/A{ro) (2.25) 

where f *^'^^(r; rg) is the substrate potential at a distance r of a sphere and a cyhnder with 
radius tq, respectively, and with 



^2.26) 



where Vexc is the excluded volume due to the repulsive parts of the interaction potentials 
(see Fig. 1). (For additional information see Ref. |l60|.) 



C. Interaction between the interfaces 



For the fiat substrate, within the sharp-kink approximation, the correction term 



uj(p){l) = is given by 



(POO POO \ 

Pi / dz &\z) -p^ dz v^'\z) , l>d^ (2.27) 
Jl-d^ Jl / 



a b 



— + — + ••• , i:^ dyj. 



Approximately the substrate potential may be expressed in terms of (Eq. (2.5)): 

v^p)(^z) = dz' J d\i <\>l^j (^^rj + j . (2.28) 

In contrast to t(p)(2;) the function diverges for z — * 0. However, according to 

Eq. (2.6), in Eq. (2.27) one has / > d^. 

In the case of a planar substrate and for the above parametrization of the interaction 
in terms of Lennard- Jones potentials the order of the wetting transition and the wetting 
temperature are determined by e, a, e^//e, cr^//a, and p^cx^j. In case of a first-order 
wetting transition a;(/) has two minima separated by a free energy barrier. One of these 
minima is localized at a finite value /q and represents the global minimum for temperatures 
T < T^u, while the second minimum is given by / = 00 where LoiV) = 0. At T = T^, the 
minimum at Z = 00 changes from a local to a global minimum, so that the film thickness 
on the planar substrate jumps from a microscopic value to infinity. In case of a second 
order wetting transition, ioiV) has only one minimum, which is shifted continuously from 
a small value /i to / = cx) as the temperature is raised from T < towards T = T^, 

m. 
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For the curved geometries, the systematic analysis of the sharp-kink approximation to 
the density functional in Eq. (2.1) leads to the result that u;{l) has essentially the same 
structure as for the planar system {I > d^): 

J^\l)^^p{^pi dr t{r;n)-p^j^ dr (^^^ v{r;ro)^ (2.29) 

with r = 1 for the cylindrical and r = 2 for the spherical system. t(r; ri) is defined 
in terms of w{r), and thus w{r), according to Eq. (2.18) and Eq. (2.19), respectively; 
v{r;ro) is the substrate potential of a sphere and a cylinder of radius Tq, respectively. It 
should be emphasized that Eq. (2.29) is vahd for any functional form of w{r) and v{r), 
provided only that they decay rapidly enough in order to guarantee the convergence of 
the integrals; Eq. (2.29) also docs not depend on the explicit form of Eq. (2.2). 

However, due to the unlimited increase of A{ro + 1) / A{ro) (see Eq. (2.13)) in Eq. (2.10) 
even at coexistence, i. e., Aflj, = 0, the film thickness / cannot take on a macroscopic 
value on cylinders or spheres, so that first order wetting transitions are reduced to a jump 
from a small to a large but finite value oil. In order to be able to discuss these transitions 
quantitatively the general expression for the effective interface potential (Eq. (2.29)) must 
be evaluated for specific interaction potentials w{r) and v{r). According to Appendices 
A and B the model potentials in Eqs. (2.3)-(2.5) lead to (see Eq. (A13)) 

u;^^\l)^u;^i{l)+ 411(1) (2.30) 

where 



rep 



UJ 



is) 
attr 



+^App,ea'^ - Z - ,f + ^° 1 ^ + I ( - ) 1, (2.31) 
90 V(^-^o)^ {h + ro)y - ll.il' V ^ 




ApA-'^^- + On-l I. (2.32) 



2^ 



(see Eqs. (B2), (B20), and (B22)) 

= ^Ap(p,ea- - p.6.;a- )^ ^; 2; (^)^ 

77r\„_^ 12^-0 ^ /II /ro\2\ ^ 99 ^ /13 U_^_fro\^\\d, 




(2.33) 
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and 




(2.34) 



Here we have assumed that the radius of the substrate is large compared with the diameter 
of the fluid particles. Equations (2.32) and (2.34) can be compared with other approaches 
p3| , |3^ , in particular with the dispersion theory developed in Ref. [34(c)] which considers 



a system composed of three concentrically arranged media characterized by different di- 
electric constants e^. The summation of the electromagnetic surface normal modes yields 
the dispersion function. Its expansion in powers of the dielectric reflection coefficient 
= results in a factorization of the dispersion energies into a frequency-dependent 
and a geometrical factor. The geometrical factors (Eqs. (6.5) and (6.9) in Ref. [34(c)]) are 
in full agreement with our results for the attractive contributions (Eqs. (2.32) and (2.34)). 
It is reassuring to see that despite of the differences between these two approaches as dis- 
cussed in the Introduction both yield the same functional dependence on the film thickness 
for large /. Furthermore one should keep in mind that the density functional approach 
yields a well-defined expression of uj{l) even for small values of /, whereas the dispersion 
theory leads to diverging energies in this limit. In addition the density functional can 
capture subdominant contributions to u!{l) which arise both from the repulsive parts of 
the interaction potentials and from detailed structures of the substrate determining the 
wetting behavior at coexistence. This will be discussed in the following section. 

For the model defined by Eqs. (2.4) and (2.5) the effective interface potential of the 
curved substrates, uj^'^\l, tq) = u'^f-lri^^ ''"o) tq), can be related to that of the planar 

substrate as follows: 

ro) = ^^Im^i (y) + a;S(0'5S (y) (2-35) 
where the scaling functions 

^ ^ ^ " -Sl(/) ' ''-H J - ^'-^^^ 

turn out to depend indeed only on the ratio x = tq/I so that 



S^i{r,/l) + ^Sill{T,/l) 
o;M(/,ro) =a;(^nO ^i^r? ■ (2-37) 



ID 

attry-) 
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Since ujfflj-.{l) = a/P and ujrep{l) = c/l^ with the Hamaker constant a = 7rAp(p^e^/(T^j — 
piea^)/3 and with c = -k /\p{pita^'^ — p^eu,/cr^^)/90, one has a;rep(/)/'^ittr(0 — c/{al^) and 
thus due to UJ^\l) = uj^lil) + Uej(0 

a;(^)(/,ro) = u;i^\l) {sj^l + (^^ (^) - 5^;^ (^))) + O(r-). (2.38) 

Therefore for Z » cr the effective interface potential attains the scahng form 

J^\l,ro) = J^\l)Sr (y) , / » a, (2.39) 

where St{x) = S^l^{ro/l). These scahng functions exhibit the following limiting behavior: 

Srix^oo)^l + ^ + 0{x-^) (2.40) 

and 

Sr{x ^ 0) = arx + (3rX^ + O {x^) (2.41) 

where ac — — — Q;^ = y, and /3s = —16. The full scaling function S-,-{x) is 

shown in Fig. 3(a). It exhibits a non-monotonic behavior. 



III. Wetting behavior 



A. Film thickness at coexistence and above the wetting transition temperature 

Wc consider a curved substrate- fluid system at a temperature T > Tw\to) above its 
(quasi-) wetting transition temperature Tw\rQ) and at liquid- vapor coexistence. (Typi- 
cally one has Tw\ro) > Tw\oo) = Tw \ see below). The curvature prevents the build-up 
of a macroscopically thick wetting film as it would form on the corresponding planar sub- 
strate. If this film thickness is sufficiently large it can be determined from the leading 
terms in the effective interface potential (a > 0): 

0.(/) c (^) + (l + + aUro). (3.1) 

Due to T > Tw^ and p — jj,o{T) = 0~ the equilibrium thickness corresponding to Eq. (3.1) 
is the maximum thickness the film can attain and it is given by 

^^L = ^T^o (3.2) 
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where x*{k) is the solution of the imphcit equation 

1/3 



X ^ K l^-j Sr{x) (3.3) 

with (see Eq. (2.36) and Fig. 3(b)) 

sAx) = (1 + x-^)(--^)/^ (sAx) + |5;(x)) . (3.4) 



This shows that the hmiting thickness of the wetting film above is governed by the 
dimensionless parameter 



(p)X V3 



^=\'^\ . (3.5) 

Due to Sr{x -^oo) = 1+ {l-l)x-^+ (^1-1 + ^^ x-^ + 0{x-^) one finds 

xliK ^oo)^(^^Y\+^-^ + {k-') . (3.6) 

1 /3 

In the limit x ^ one has Si{x ^ 0) = (^) x~^/^ (l + ^x — ^x"^ + 0{x^)) and 
s^{x ^ 0) = (1 + fx - fx^ + 0{x^)) so that 



<..('.-0)=fS)''Vn"Vo(.»/=). (3.8) 



and 



2/ V2/ 

This implies that for a given system, i. e., and a fixed and for large Tq the maximum 
thickness of the wetting film is given by 

1/3 / \ 



Apro 



(see also, c. f., the first paragraph of Subsect. III. B) whereas for a small liquid- vapor 
surface tension a^^ and a large Hamaker constant a the maximum thicknesses are (see 
Eq. (3.7)) 
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and 



^^^^-1 VT^) ' ' (3.10) 



For large k Imax is larger for a cylinder than for a sphere by a constant factor 2^^^ ~ 1.26 
(Eq. (3.6)) whereas for small n this ratio diverges weakly as (97r)^/^2^-'^'^/^*^/s:~^/^° (Eqs. (3.7) 
and (3.8)). The full dependence of Imlx/ro = l/x* on k is shown in Fig. 4. 

In principle k can be any positive number. Typically, at low temperatures, al^g ranges 
between 10"^ J/m^ and 10~^ J/m^ and a between 10"^°J and lO^^^J so that 
K = kq ■ (ro/A)^/^ with kq = 0.05 . . . 0.20. With kq = 0.05 according to Fig. 4 this implies 
that for To = lOOA Imax — 130A and Imax — 90A whereas for tq = l/xm Imax — 550A 
and Imax — 430A. Due to finite resolutions for pressure and temperature equilibrium 
film thicknesses on a planar substrate can be followed up to at most a few hundred A, 
say 500A. In this case and for the above quoted values of a^^g and a tq must be smaller 
than 50/im for a cylinder and lOO/im for a sphere in order to cause a deviation from the 
corresponding planar film thickness of more than 5%. 

If for a given system, i. e. tq fixed, the temperature is raised towards Tc, cr^^^^ vanishes as 
aoP", V ~ 0.632, where t = {Tc-T)/Tc 0, and a ~ Ap as aot^, /3 ~ 0.328, so that k goes 

\ 1/3 _ 

to zero as ( t(2i^-/3)/3 ^ p. 312 ^-^^ Eqs. (3.7) and (3.8) apply. From these one infers 

that the maximum film thicknesses diverge as Imax — (97rrgao/2cro)^^^^ '■^'^^^''^^ ~ ^-0.234 
and l^mlx ^ (8ri]ao/ao)^/^t~(2''-^)/^ ~ respectively. At the critical fluctuations 

in the bulk lead to a power-law decay of the density profile towards the bulk density. 
As function of the temperature T this mesoscopic interfacial structure known as critical 
adsorption emerges from a film whose thickness diverges like the bulk correlation length 
^ = ^^t'^ — ^(7^^°'^'^^ where is typically of the order of 1.5A. This picture for a 
planar substrate holds also for spherical and cylindrical substrates provided ^ is large 
compared with their radii of curvature |26|. On this basis we expect that ultimately for 
f — *• both Imax and Imax should diverge as t'" . Our above results obtained from the 
effective interface potential signal such a divergence but the corresponding exponent is 
smaller. This difference must be attributed to our present sharp-kink approximation of a 
mean-field theory. Nonetheless one should keep in mind that at least within the mean-field 
theory for a planar geometry the leading behavior of the effective interface potential is 
captured correctly by the sharp- kink approximation in spite of the actual broadening of the 



emerging liquid- vapor interface close to [55|; it is reasonable to expect that this quality 



of the sharp-kink approximation carries over to the case of curved surfaces. However, as 



15 



the discussion of Eqs. (2.46) and (2.47) in Ref. [55(c)] shows even this quahty of the sharp- 
kink approximation does not suffice to describe properly the crossover between complete 
wetting and critical adsorption. The crossover condition I ~ ^ must be inserted by hand. 
Our above findings for Imlx must be judged and augmented similarly. 

The above analysis relies on a positive value of a and, because we have considered only 
the leading terms in uj{l), on the fact that Imix is large. This is fulfilled for T > Tw\ro) 
where Tw\ro) is the thin-thick transition temperature in the case of a (quasi-) first-order 
wetting transition. If the wetting transition on the planar substrate is continuous, there is 
no transition at coexistence and we consider simply T > Tw\ i. e., a > 0. The sohd- vapor 
surface tension is the minimum value of Jls(/) so that at coexistence 

A(rQ -\- Iq) 

(Ts,g = minQs(Z) = Vts{lo) < ^/,g(^o + ^o) + 'yi,s{ro) (3.12) 

because a;(oo) = so that < (see Eq. (2.10)). One should note that in Eq. (3.12) 
the inequality never turns into an equality if Tq < oo because Zo('"o < oo) < oo so that 
uj{lQ{rQ < oo)) < 0. Only in the limit tq — > cxo one recovers the standard relation that 
Csfa < cj^g + o'j^] turns into an equality at T = Tif\ Since in Eq. (3.12) the right hand 
side ai^g is enhanced by the factor A{rQ + lQ)/A{rQ) as compared with the planar case one 
can surmise that Tw\ro) is larger than Tw^. However, this argument must be checked 
more carefully, because at T^\ro) the right hand side of Eq. (3.12) is not attained. This 
will be done in Subsec. III. C. 



B. Complete wetting off coexistence and above the wetting transition temperature 

According to Eq. (3.9) the thickness Imlx of a wetting film on a curved substrate at 
liquid- vapor coexistence A/i = is the same as if the identical system in planar geometry 

is taken off from liquid- vapor coexistence with a Afj,eff{Afj, = 0) = ':p^Ap- ^^^^ subsec- 
tion we want to discuss how a wetting film on a curved substrate attains this maximum 
thickness Imlx when one approaches liquid- vapor coexistence, i. e. A/i — 0, above Tw\tq). 
In this case we have (see Eqs. (2.10), (2.12), (2.13), and (3.1)), with a > 0, 



a(/»ci.)^^^. (y) 



1 + - 1 a 
To 



+ 



so that 
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where y*{K, A) is the solution of the imphcit equation 

/r\ 1/3 

with 

Ur{y, A) = f 1 + - fl + V sr{y). (3.16) 



T \ y 

The dimensionless parameter A measures the distance from hquid- vapor coexistence, 

A=^A/., (3.17) 
^1,9 

so that y;{K, A = 0) = x;{k) = ro/lmax{K). With Uriy ^ 00, A) = (l + J)^^^ + y-^{T + 



A - |)/4(1 + A/t)2/3 + 0(y-2) one has 

y;{K - 00, A) = (^) « (1 + ^) + + O {k-') . (3.18) 



In the hmit k one finds Ui{y ^ 0) = {^f^^ (^1 + + O ((f)^)) and 

-2(y - 0) = (A)V3 1(1+^-^1 + ((|)^)j so that 



yr(.^0,A)=(^j +^(^) +0(.^/^) (3.19) 



3A\'/' 1 + 3A /k'xY^' 



and 

^T6"J ^ 3A [1^) 
In the case k » 1 Eq.(3.16) imphes 



y;(«-0,A)= ^ +0{k'/'). (3.20) 



/o(ro, A/.) = 3 , fi; » 1, (3.21) 

(1 + 7)^ 



whereas for k <^ 1 one has 



and 
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Zf)(ro,A/.)=(^y^', (3.23) 

Equations (3.16), (3.22), and (3.23) show that the hmits A — > and k — > cannot be 
interchanged. 

Since for a small undersaturation A/i the thickness of the complete wetting film on the 
corresponding planar substrate is (T > Tw^) 

Zo(ro = oo. A/. ^ 0) = l^f\A^, ^ 0) = 1^-^) , (3-24) 

the film thickness on the corresponding curved substrate can be expressed in terms of an 
effective A//e//(A//; tq, r): 

/o(ro. A/. ^ 0) = /?^(A/ie//), (3.25a) 

A/ieff ^ Aii + Sii{k,X). (3.256) 

In the following we determine A/Xg// on the basis of Eq. (3.13) which captures the leading 
behavior of Qs{l) for large /. (One should keep in mind that prewetting phenomena depend 
also on the behavior of Qsil^dw) so that one cannot expect that the (quasi-) prewetting 
line on curved substrates is simply the prewetting line of the planar substrate, A/i,p^(T), 
shifted upwards by Afx^ff- Thus Eq. (3.25a) holds only above the prewetting hue where 
lo is sufficiently large. Therefore we concur with Upton et al. (Sect. XI in Ref. [26(b)] 
in that the effects of curvature cannot be fully subsumed into an effective A/Zg//.) With 
Eqs. (3.14) and (3.17) one has 



Ae,,(.,A) = ^^^^ = 2(^)\ (3.26) 



In the limit k ^ 1 Eq. (3.18) apphes so that Aji^ff is given by 

(p) 

A//e//(AAt;ro,r) = A/x + ^^, « > 1. (3.27) 

Thus in the limit k ^ 1 the curved substrates resemble the same film thickness as the 
planar substrate if the liquid- vapor coexistence curve iJ>o{T) is shifted upwards by 



^/^^^ = ^- (3.28) 



ZI2 

Apro ~ ApaV2 

This is in accordance with our previous finding Ajjif.ff{AjjL = 0) = Sjj!'^ (see Eq. (3.9)). 
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From Eq. (3.27) one finds that for spfieres (r = 2) tfic effect of tlie curvature is twice as 
large as for cylinders (r = 1). It should be emphasized that the shift is independent 
of the Hamaker constant a and depends on surface properties only via the planar liquid- 
vapor surface tension. 

Equation (3.27) holds only in the limit = (aJ;VoVa)V3 » 1, i. e., for To — > oo and T 
fixed. In the opposite case, i. e., for T ^ Tc and ro fixed one has k ~ _> g so that 

the limit k <^ 1 must be considered (see Eqs. (3.19) and (3.20)). These latter equations 
are valid provided that (i) y* <^ 1 and (ii) (l + l/y*)A/r ^ 1 (compare Eq. (3.16)). 
Condition (i) is fulfilled if k < (97r/2A)^/^ for a cylinder and if k -C for a 

sphere, respectively. Condition (ii) is always fulfilled for A > r; for A < r it requires 
y* < A/(t-A), i. e., K < (97r/2)i/=^A^/V(l-A)5/3 for a cylinder and k < 2^/^\^/y{2-Xy 
for a sphere, respectively. In this regime, denoted as I (see Figs. 6(a) and (b)), Ag// is 
given by 

(2A\ 
— j regime I, (3.29) 



and 



AX V2 



Xl;jf - ( 7 ) regime I. (3.30) 



On the other hand according to Eq. (3.18) the limit y* ^ 1 applies ii k,^ 2^^^/ (r + A)^/^ 
(regime II); in this case Eq. (3.26) yields 

Ag^j ~ r + A, regime II, (3.31) 

which is equivalent to Eq. (3.27). In the intermediate regime III, i. e., 1 ^ y* ^ X/ (r — A), 
the asymptotic formulae in Eqs. (3.19) and (3.20) do not hold. Instead of expanding y* 
in powers of k, in this case one has to evaluate Eqs. (3.15) and (3.16) for small values of 
A first and then to consider the limit k ^ 0: 

y; {k,X^O)c^k [^Y' (l + ^ (l + ^)) 'r{y). (3.32) 
With Sr{y) taken in the limit y ^ one finds in this regime III 



and 




fl-liM I 1 + 7 I (i^) +2 + 0{^'^))+0{>?)) (3.33) 



- iT + To ((i)"'^ ^^"'^ + ° (^^)) • P-^) 
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respectively, so that 

2 V^U. 3. ( f2K'\-"' 



and 



regime III, (3.35) 



1/5 



3 , f fK\-^l^ 



regime III. (3.36) 



Figures 7(a) and (b) visualize this crossover behavior of X^jf as function of k for A fixed 
in the case of a cylinder and a sphere, respectively. 



C. First-order wetting transitions 



The results of the previous subsection allow us to discuss the effect of curvature on 
a wetting transition which is first order on a planar substrate. For suitably chosen in- 
teraction potentials Fig. 8 shows such a phase diagram and the temperature dependence 
of the corresponding film thickness at coexistence for a cylinder and a sphere, respec- 
tively, as compared to those of the corresponding planar substrate. In accordance with 
Refs. pB| , |5^ we find for the planar substrate that the prewetting line /ip^(T) joins the 
gas-liquid coexistence curve fioiT) tangentially, 

/x(^j(T^T(^))-/io(T) = A;,|r-r(^)r(^), ^(C) = (3.37) 

if for large distances the pair interaction potentials in the system decay like r"^^"*"''^ ||63|| . 
Upon crossing the prewetting line the wetting film thickness undergoes a thin-thick tran- 
sition where the corresponding jump Al{T) diverges at the first-order transition temper- 
ature Tw^ as 

Al (T ^ Tif)) =ki\T- T^r^^'^)/'^. (3.38) 



In Subsects. III. A and III.B we focused on the region T > of the phase diagram 
where the behavior of the film thickness is independent of the order of the wetting tran- 
sition. Since, however, in the region close to Tw^ the order of the transition plays an 
essential role, the present and the following subsection are dedicated separately to first- 
and second-order transitions, respectively. A first-order wetting transition is accompanied 
by a prewetting line which implies the occurrence of discontinuities in the film thickness 
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also off coexistence. In tfie region between tlie prewetting line and the coexistence curve, 
i. e., T > and A/x < Anp2{T) (see Subsect. III.B), the film thickness on the curved 
substrates and on the corresponding planar substrate are related to each other via A/^g//- 
In the limit of large radii, k ^ 1, this relationship is given by the shift AyUe// = A// + S/i. 
Thus one is inclined to surmise that for k ^ 1 also the prewetting line of the curved 
substrates can be obtained by shifting the prewetting line of the corresponding planar 
substrate upwards in the phase diagram by the amount S/i. In this case the knowledge 
of the location of the prewetting lines relative to each other yields information about the 
shift of the transition temperature, Tw\ro) — Tw \ and the jump in the film thickness 
AZ^'^) on the curved substrates. Since Eq. (3.37) holds especially for T — Tw\k oo), 
one has with ij,^{Tw\k oo) ^ Tw^) — ^q{Tw\h oo)) = 

/ 1 \ 1/7(0 

tI^\k ^ oo) - Tif) = [j-^^^t^j (3.39) 
so that for the jump in thickness at coexistence on the curved substrate one has 




(3.40) 



Thus C = 3 implies Tir^(«; ^ oo) - T^^ = (37rT/(4A;^))3/2(Ap/ro)3/2 and M^''\T^\k 
oo)) = fcK4V(3vrr))V3(ro/Ap)V3. 

However, the results of Subsect. III.B hold only for large / whereas the study of the 
prewetting line involves also the behavior of VLa{l) for / ~ dyj. This requires to test the above 
considerations numerically which takes into account the full function VLs{l). Therefore we 
have calculated the prewetting lines for the planar, spherical, and cylindrical geometries 
for fixed fiuid-fiuid and fiuid-substrate interaction potential parameters and for different 
radii. The analysis of the dependence of the jump in the film thickness and of the shift 
of the transition temperature yields {Tw\ro oo) — Tw^)/Tw^ ~ f^o"^^^ ^^"^ Al^'^\ro — > 

1 /3 

oo) ~ Tq which confirms the above heuristic considerations. However, the inspection 
of the prewetting lines in Fig. 8(a) demonstrates that the phase diagrams of the curved 
systems do not correspond to a simple shift of the phase diagram of the planar system. 

Since k = {9{^)^Ap/{4{p^C-^){^)^ - pi)))^^\ the limit k < 1 requires (i) a small 
radius of the substrate, (ii) large substrate potential parameters p^ewfa^j^ in comparison 
to pi, or (iii) T close to Tc. Simultaneously the parameter space of p^, e^f, and a^f is 
restricted by the requirement to have a first-order wetting transition in the temperature 
interval between the triple point temperature T3 and Tc. Within the framework of our 
numerical analysis a set of parameters which fulfills these conditions has not been found. 
Thus a numerical analysis of the case k <^ 1 would require to enlarge the parameter space 
considered here which we have not pursued. 
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D. Critical wetting transitions at coexistence 



Up to now we have considered the case that the Hamaker constant a of the planar 
system (Eq. (2.27)) is positive. This is always valid for T > Tw \ If the wetting transition 
is first order this holds even for T = Tw \ In the latter case the description of the system 
for T < Tw^ requires one to consider the full effective interface potential and not only its 
leading behavior for large film thicknesses (Eq. (2.27)). 

If, however, the substrate potential is sufficiently weak, a favorable choice of interaction 
potentials allows for the occurrence of a continuous, so-called critical wetting transition 
such that a changes sign at Tw^ and the coefficient b of the next-to-leading order term in 
ti;(/) is positive at that temperature |^ (Eq. (2.27)): 

a(T ^ T^) ^ 0, a{T = T^) = 0, b{T^) > (critical wetting). (3.41) 



For critical wetting there is no prewetting line. Recently such a continuous wetting tran- 
sition at coexistence has been observed experimentally |^ . This implies that on a curved 
substrate any wetting transitions, which are continuous in the planar geometry, are wiped 
out. Due to the effective upward shift of the bulk phase diagram caused by curvature (see 
Subsect. III.B), the wetting behavior on a curved substrate upon raising the temperature 
along the liquid- vapor coexistence curve = fioiT) — is that of the corresponding pla- 
nar substrate if in the latter case one passes the critical wetting transition at Tw^ with 
a constant undersaturation 6jj,'^ parallel and below ij.o{T) (Eq. (3.28)). This results in 
a steep but smooth increase of the thickness of the wetting film on a curved substrate 
upon passing Tw^ at Ho{T) and a levelling off at a finite value for T > T^^ (see Fig. 9) 
for which the previous considerations (Subsects. III.B and C) hold. Since, however, these 
previous results have been derived for the case a > 0, this analysis must be refined in 
order to cover the case that a changes sign at T^. 

To this end we first consider the circumstances under which critical wetting transitions 
can occur. As mentioned above, the substrate potential must be sufficiently weak so that 
a(T) < for T close to the triple point. However, for pure Lennard- Jones interactions (see 
Eqs. (2.4) and (2.5)), it is well known that the excluded volume effect near the substrate 
(i. e., dw > 0) leads to b^T^j) < and thus induces a first-order wetting transition pil|J55| . 



One can fulfill the necessary condition b{T^) > (Eq. (3.45)) by coating the substrate 
with an appropriate thin film consisting of a different material which leaves the leading 
asymptotic behavior of the substrate potential and thus a unchanged but introduces a 
suitable subdominant contribution which can overcompensate the effect of the excluded 
volume within the effective interface potential. 

Therefore we consider a planar substrate filling the half-space z < and consisting of 
two types of particles which are located in the region z < —6 and in the film —6 < 2; < 
and which interact with the fiuid particles according to Eq. (2.5) with interaction potential 
parameters {o'^j,^^^}} Wwj^^wl}^ respectively. With the simple assumption that 
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the two species exhibit the same number density pw this leads to a substrate potential 
(Eq. (2.28)) V^P\z) = p^v^p){z) with 



Us U4 

— + — + ... 

z^ . 



(3.42) 



where 



"3-3 ^yjf [^wf J 5 



(l)^Jl)^6_ (0). (0)x 

wf) '^wfK'^wf) 



(3.43) 



(This expression for does not yet include contributions due to a discrete arrangement 
of the substrate atoms on a lattice; in this respect see Ref. [0.) Thus with (see Eqs. (2.4) 
and (2.15)) 



+ 0{z- 



(3.44) 



and 



27r 



-ecr 



(3.45) 



one has 



a = -Ap {u^pw - hpi) 
b = ^Ap {u^pw - Sdy^tspi) 



(3.46) 
(3.47) 



Thus 6 > is realized for a sufficient large value of M4 which is, e. g., the case if the 
coating of the substrate consists of particles whose interaction with the fluid particles is 
stronger than that of the particles in the main part of the substrate beneath. 

Analogously, for curved substrates we consider coated spheres and cylinders which 
consist of particles characterized by interaction potential parameters {e^]-,cr^°j} and 

{e2^|,o"^^]} for < r < ro — 5 and tq — 6 < r < r^, respectively. Based on the results 
of Appendix C Eq. (2.35) can be generalized in order to be able to describe continuous 
wetting transitions: 



-;^^'(t)4'^'(t) 



(3.48) 



where the scaling function 7lr{x) (see Fig. 10 and Eqs. (C4)-(C7)) is given by (x = tq/I) 



TZr{x) = lim 

i— »oo 



73 

b 



(3.49) 
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For large 7^^(a;) takes on the form 



Urix ^ oo) = 1 + ^ + 0{x-'^) (3.50) 

whereas for x ^ one has 

n^^x ^ 0) ^ — {1 - 3x + 0{x'^)) (3.51) 

and 

7^,(x^0) = 8(1- 3x + 0(^2)). (3.52) 
The minimization of the corresponding effective interface potential 

a.if) = ^5. (^) + (^) + (l + ;^) ' + (3-53) 

yields for the equilibrium film thickness the implicit equation (note that here «; < is 
possible) 



with 



and 



r^(a;, 7) = (1 + 7x5^(0;)) ^^^s^(a;) (3.55) 



The functions Sr{x^ and gr(a^) depend only on the geometrical properties of the substrate 
whereas the fluid-fluid and fluid-substrate interaction potential parameters enter in form 

of the dimensionless numbers k — {r^aflla^ and 7 = 36/ (2aro). For the explicit forms 

of the functions q-iix) see Eqs. (C6) and (C7). The knowledge of the behavior of qri^x) 
and St-[x) in the limits a: ^ and x ^ 00 allows one to solve Eq. (3.54) approximately 
for specific regions of the (k, 7) parameter space (see Fig. 11). Since, however, one has 
X = Tq/I > 0, only the regions (k < 0, 7 < 0) (for T > Tw '') and (k > 0, 7 > 0) (for 
T > Tif^) are of interest. 

With qr{x ^ 00) = 1 + \x-^ + ^x-^ + O {x'^) and s^(x ^ 00) ^ 1 + [\ - \) x'^ + 
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^ X + O {x ^) one has for I7I 3> ^: 



24 



SO that the solution 7) of Eq. (3.54) reads (note that is always positive) 

x:i\K\ 00, i7i - = Y j " 4^ + 2 - ^)) + ^ (M"') ■ 



(3.58) 



Thus in the limit of fixed but large ro and a ^ one finds for the film thickness ij^^ — ro/x* 
{h must be positive for critical wetting to occur) 

/ X 1/4 / X 1/4 /X 1/3 

where is defined by Eq. (3.28). For the corresponding planar substrate the film 

thickness along a complete wetting path A// — > approaching the transition temperature 
T = Tw^ of the second order wetting transition is given by /o(?"o = cxd, A/i, T = Tw^) = 
( A&)^^^ = li^\T = Tir\ A^) due to a(T = T^f^) = (see Eq. (2.27)). Thus in the limit 
considered in Eq. (3.58) the film thickness at coexistence on a curved substrate equals the 
film thickness of the corresponding planar system which is off coexistence at Ajj, = djj,'^ . 
This imphes that the picture, which has emerged in the previous subsection for a > 0, 
carries over to the case a — 0. 

For a; — > 00, I7I <S I the sign of 7 has to be considered. If 7 is positive, i. e., T > 
one can expand r^(x) in powers of as before, which leads to 

rr{x ^ 00, 7a: < 1) - 1 + ^^^2:-' + ^ + ^^6"^"^ ' ^ ^ ^^~^)) 

+ C(x-^7V) (3.60) 

so that one finds for the solution of Eq. (3.54) 

'r\V3 3r-4 



+00,7^0+) = fi:(^-j + 
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Hence for large Tq and b/a ^ the film thickness attains its asymptotic value given by 
Eq. (3.9) as 



3h 

T>Ti^\ ro»->0. (3.62) 
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For 7 — s> , K — s> — oo, i. e., for T < Ti^\ ro ^ oo, and a fixed, the film thickness is 
determined by the equation 

+ = ^K^. (3.63) 

Since 0:^+70;^ has zeros at a; = and x = — 1/7 and a maximum at x = —3/ (47), Eq. (3.63) 
implies x* > — 1/7 = Tq/I^q^ and thus Z^"* > With the ansatz x* — — 1/7 + 5x 

Eq. (3.63) leads to 



X 

7 



-i(l + ^7V + 0((7«:)')), 7'^«1, (3.64) 



so that one has for the film thickness 



lt\T)^l^^ (T, A/. = 0)11 



16 roo^ 

1/3 



T<Tl^\ ^o» . (3.65) 

So far we have discussed the limit x — > 00. For x ^ one has Sc{x — > 0) = 
(^^) (1 + 1^ - i^' + ^(^'))' ^^(^ ^ 0) = ^ (1 + fa: - + 0{x^)) and g,(x ^ 
0) = A(i + |a;2 + 0(a,3)^^ ^^(3, ^ 0) = i (1 + 2,t2 + C(a;3)). Thus r,(a; ^ 0,7) 

is given by r^x 0,7) = (^)^^^ (l + §7)'^^^ ^ (l + fa^) and rs(x ^ 0,7) = 
1(1 + 7)"^/^ ^ (1 + l^) which leads to 



for a cylinder {t — 1) and 



for a sphere (r = 2), respectively. Obviously for a < there exists a physically meaningful 
solution only if 7 < — The film thicknesses are thus given by 



Qtt arn b 



1/4 / >^ 1/2 

lal 



and 
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;(«) 



1 + 



36 
2aro 



1/5 



1/2 



(3.69) 



Upon crossing Tw^ the parameter 7 jumps from —00 at Tw^ — to +00 at Tw^ + 
because at T = Tw^ the Hamaker constant a changes sign. At Tw^ x* depends only on 
the parameter = 2rQal^g /{3b) which is positive both above and below Tw^: 



:(H«i.l,Hoo)^(^)"\(^)"%o(( 



.3 \ 1/2 



X. 



1^ 

— <1, 

7 



,3\ 3/4^ 



(3.70) 



and 



/^3\l/5 / 

x:(i«i«i,i7i-oo) = (^-j +(^ 



K 



7 



«1, 



^3\ 2/5 



so that one obtains for the film thicknesses at Tw^ the results 



^3\ 3/5^ 



lie) _ 
In — 



(drcbrQ 



1/4 



36 



1/3 



2a 



(p) 
1,9 



(3.71) 



(3.72) 



and 



lis) 



12ro6 

^9 



1/5 



36 



1/3 



. 2(T; 



',9 



(3.73) 



In the limit 7 0^, i. c., tq fixed but small and b/a — > 0, one recovers the first-order 
results given by Eqs. (3.10) and (3.11), respectively: 



/f)(|.|«l,7-0-)=ro (^) 

r>Tir), 0<^«^ 

a 3 



1/4 



1 + ^7 + (7^) + O i^y^) 



(3.74) 



and 



/(^)(l'^l«l,7-0+) = ro(-) 



2\ 



3/5 



K, ) 



(l + ^7 + 0(f))+0(«^^/^) 



T > T(P) < - < — 
' a 3 ■ 



(3.75) 
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For those regions of the (k, 7) parameter space which are not covered by the above con- 
siderations of the hmits x — and x ^ 00 the solution of Eq. (3.54) requires a numerical 
treatment of Eqs. (3.54) - (3.65) using the full expressions for the functions Qrix) and 

Sr{x). 



IV. Finite-size effects 



In Sect. Ill we have shown that the positive curvature of a substrate prevents the forma- 
tion of infinitely thick wetting films. At coexistence a critical wetting transition is smeared 
out whereas a first-order wetting transition is reduced to a discontinuos transition between 
a thin and a thick, but finite, film. However, the latter nonanalytic behavior as well as 
the discontinuities associated with the corresponding prewetting line are artefacts of the 
mean-field theory we have used. Since the surface area of spheres is quasi-zero dimensional 
and that of cylinders quasi-onedimensional, fluctuations erase any phase transition and 
lead to a rounding of the above mentioned discontinuities. With increasing radii of the 
spheres and cylinders this rounding sharpens and a first-order transition emerges gradu- 
ally. In this sense our mean-field results become increasingly more accurate as the size of 
the spherical and cylindrical substrates becomes larger. 



Based on the results of Privman and Fisher for the influence of finite-size effects 
on first-order transitions we are able to estimate the rounding of the first-order wetting 
transitions described above (see also Refs. [^] and |]2B|). In Ref. |]5^ the scaling functions 
for the singular part of the free energy and the magnetization in an external field of a finite 
Ising spin lattice, i. e., a stripe and a block geometry, are derived. The shape of these 
functions allows one to determine the region over which the phase transition is smeared 
out. In order to translate the results of Ref. onto the case of the prewetting transition 
we analyze the effective interface potential ^2^(0 given by Eq. (2.10). In the vicinity of 
the prewetting line fipre{T) fis(/) exhibits two local minima, l^{T,fi) and ly{T,fi) with 
the corresponding values fis(/<(T, /x)) =: flf{T,fj,) and Qg {i>(T, fi)) =: Qf{T,fi) of the 
effective interface potential. One of the two thicknesses /< (T, /i) and /> (T, fi) corresponds 
to the global minimum which represents the equilibrium film thickness l^^\T, fi) at the 
point (T, /i) of the phase diagram. On the prewetting line one has Qf{T, Hpre{T)) = 
^1^{T, jj,pre{T))) so that the thin film {1^^\t, jj,pre{T))) and the thick film {1^\t, jj,pre{T))) 
coexist. We now consider the difference AQs(T,fi) = f2^(T, /x) — Qf(T,^) between the 
two minimum values for a given temperature T as a function of the chemical potential /x. 
According to Ref. |^| the width Sfi{T) of the rounding region for a spherical substrate is 



given implicitly by the inequality 

\Ans{T,^i = fipre{T) + 6fi{T))\rl ^ ksT. (4.1) 

On a cylinder along its axis the liquidlike film is expected to break up into domains with 
small (/<■*) and large (/>"*) film thickness, respectively. The corresponding kinks at the 
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domain boundaries are accompanied by an additional free energy per unit length, i. e., 
a line tension S;. Following Ref. we estimate Hi by the product S; ~ Ala 



(p) 



Ap) 



l>(yf'g- Reference provides a more detailed picture of the line tension 



(^> - ^<) ^l,g 

associated with the coexistence of thin and thick films on a planar substrate along the 
prewetting line. We note that from these analyses one expects that in systems governed 
by nonretarded van der Waals forces E; diverges logarithmically upon approaching coex- 
istence at Tyj. Since for large tq one has A/ ~ |T — T^|~2 (Eq. (3.38)), our ansatz for 
overestimates the actual value. Since we aim only for a rough estimate we refrain from 
a more detailed analysis along these lines. (For a discussion of finite-size effects on the 
line tension see the work by Indekeu and Dobbs in Ref. |6^.) The rounding region of the 
prewetting phase transition on a cylindrical substrate is then given by |66 



|Af],(r,/x = /ip,e(T)+5/i(r))|roeexp(roSKT,/i = /ip,e(T))/fcBT) ^ keT, (4.2) 



where ^ is the bulk correlation length. The width Sfi{T) as determined by Eqs. (4.1) and 
(4.2) can be translated into a temperature interval 6T{T) by taking into account the slope 

m^re{T) = d (/ipre(T)/e) /O {ksT /e) of the prewetting line: 



A;B5r|^|m(;)(T)r'|5M2^)|. 



(4.3) 



Now consider AQg at a point (T, fipj.e{T) + 6fi{T)) on the boundary of the rounding region 
given by Eqs. (4.1) and (4.2), respectively. The calculation of AQs{T, ^pre(T) + 6fi{T)) 
requires the evaluation of the effective interface potential for the two minima at /> and /<. 
Whereas in the asymptotic form (Eq. (3.13)) of Eq. (2.10) yields sufficient information 
about the second minimum Z>, this is less transparent for /<(T, fi) and ^f{T, n). However, 
because Qf (T, /i) varies only weakly upon a deviation 6fi from the prewetting line ^ = 
fipre{T) [H, one has Qf {T,fipre{T) + 6fi) ^ Qf {T,fipre{T)) = Q>{T,fipre{T)). Thus one 
obtains 



{T,fXpreiT)+Sfx)\ ~ \Sfx\ 



dn>{T,fi) 



(4.4) 



fJ.=fJ.pre(T) 



With Eqs. (4.1) - (4.3) this implies 
\kBoT\ ^ — 2- 



\m. 



is) I 
pre I 



dn>{T,fi) 



(4.5) 



A'=Alpre (T) 



for a sphere and 
\kB5T\ 



ksT 



\m. 



pre \ 



exp (-roS^/ZcfiT) 



(4.6) 



/^=Mpre{T') 



for a cylinder, respectively. The factor dQf /dfi consists of four contributions: 
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mi 

dfj, 



+ 



i->,Pl,Pg 



+ 



mi 

dl> 

mi 

dpi 



P,l>,Pg 



dl> 
dfj, 



m; 



P,l>,Pl 



(4.7) 



Since the variation of = /iq) - pg{iJ,)) / pg{p = po) and {pi{p = po) - pi{p)) / pi{p = po) 
are both of the order 10~^ along the whole prewetting line, the last two terms of Eq. (4.7) 
can be neglected. With Eq. (3.13) for the effective interface potential fls{l) one obtains 



mi 

dp 




T+l 



r-1 



1 + A 

T 



1 + 1 

y 



with y — Tq/I^. The second minimum Z> is given as the solution of the implicit equation 

(3.15). Depending on the value of the dimensionless parameters k — (j'ocri^g /a^ and 

A = (36/(2aro)) there are three different regimes in each of which the implicit equation 
(3.15) provides an approximate solution /> (see Sect. III.B and Fig. 6). Accordingly, this 
leads to three different limiting cases of the estimates (4.5) and (4.6) for the rounded 
regions. 

1. In region I one has k ^ 0, X fixed, and y given by Eqs. (3.19) and (3.20) which 



leads to 



(9/1 



9/1 



-1 - 

— I ^^^^^ ^ for a sphere and 
cylinder, respectively. Thus with A/Zp^e — Po ~ Ppre one finds for the width of the 
rounding region \kBST\: 



\kBST\ 



3 keT 
2 



Tflpre I 



r^aApJ 



1/2 



(4.9) 



for a sphere and 

\kBST\ 



X cxp 




(4.10) 



for a cylinder. 

2. For film thicknesses /> small compared with the radius Tq (region II) one 
has K ^ oo, A fixed, and y is given by Eq. (3.32). With 



roAp I T roa 



1/3 



1/3 



-1 



1 I rpApA^i 



one finds 



30 



\kBm 

for a sphere and 

\kBm 



2a 



(p) 



ro 



Ap\m 



pre I ' 



2a 



1/3 



pre 



(4.11) 




(4.12) 



for a cylinder, respectively. 

3. In region III one has A/k — > and y given by Eqs. (3.33) and (3.34), so that one 
finds 



\kB6T\ ^ 
for a sphere and 

\kBm - 



ksT 1 / (T'l 



(p) \ 



,9 



1/5 



1 + — 



3 / Sarf, 



\ 



10 \ Jp) 

^1,9 



6 



ApAfi. 



pre 



(4.13) 



fc^T 1/2 fx,,. 



X exp 



1/4 




(4.14) 



for a cylinder, respectively. 

From Eqs. (4. 9)- (4. 14) one infers that for a sphere the width of the rounded region 
shrinks for tq — > oo always algebraically but exponentially for a cylinder. In the case 
of large radii and small film thickness (region II) our result for ST for a sphere agrees 
with Eq. (5.5) in Ref. p5[. As already stated in this reference it is difficult to obtain 
a general expression for ppre{T) and mprL However, our microscopic model allows us to 
obtain quantitative estimates for ppre{T) and m^re for specific interaction potentials of the 
system. From the first-order wetting phase diagram of a cylinder and a sphere of radius 
ro = lOOo" (see Fig. 8) one can infer A/Zp^e directly as a function of the temperature. If 
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one chooses /c^T = 0.95e, one has ksTw^ ~ 0.785e, AfipJe — O.Ole, {mpXl ^ — 1.65, and 
Afi^pl ~ 0.026e, \m^^e\~^ ^ 1-45, respectively. Thus with p^a^ (fc^T = 0.95e) ~ 0.076, 
PiO'^ {ksT — 0.95e) ~ 0.551, as well as the substrate parameters e^/ = 0.625e, a^f — 
1.35(7, and p^a^ = 1.0, Eq. (4.11) yields 



\ST\ 



J- 1 



2.3 • 10-^ (sphere, region II). (4.15) 



Equation (4.12) for a cyhndrical substrate involves also the bulk correlation length ^ de- 
fined as = I / (frr'^Qir)/ J (PrQ{r) where Q{\r — r'\) =< /3(r)p(r') > — < /5(r) > 
< p(r') > is the two-point correlation function of the fluctuating number density p(r). 
Within our model this leads to = ^\w2\p'^k,t with the isothermal compressibility 

kt = {l)^ + Wo j and Wi = f d'^r rhv{r). For the model potential of the fluid- 
fluid interaction given by Eqs. (2.3) and (2.4) one obtains ^{T = 0.95e/kB) — l.Tcr which 
implies 

\ST\ 



T 



10~^*" (cylinder, region 11). (4.16) 



We conclude that for ro = lOOcx the rounding of the prewetting transition on a sphere 
is at the lower end of the experimentally accessible temperature resolution whereas for a 
cylinder the rounding cannot be detected. 



V. Summary 



We have studied wetting phenomena on cylindrical and spherical substrates (Fig. 1) for 
simple fluids whose particles are governed by dispersion forces and are exposed to long- 
ranged substrate potentials. Our approach is based on a microscopic density functional 
theory. Using a sharp-kink approximation for the density proflle (Sect. 11 A), we have 
determined the effective interface potential Q-sQ) (Eq. (2.10)) for the emerging hquidhke 
film of thickness I. From a detailed analytical as well as numerical discussion of fis(Z) for 
a cylinder {t — 1) and a sphere (r = 2) of radius tq we have obtained the following main 
results: 

1. As a contribution to VLs{l) we have determined the liquid- vapor surface tension 
a\^g (/i) of a spherical drop and a liquid thread of radius h (see Subsec. II B) . For small 
values of |- it approaches its planar value af'^ = (^'fgih — oo) nonanalytically: the 
next-to-leading order term beyond the contribution ~ ^ is proportional to ln(/i) / h'^ 
(see Eqs. (2.21) - (2.23) and Fig. 2). 
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2. The interaction contribution uj^'^\l,ro) to the effective interface potential fls{l) has 
been determined in closed form (Eq. (2.29)). It can be expressed in terms of the 
corresponding planar quantity u;^\l) (Eq. (2.27)) by scahng functions Sr{ro/l) and 
TZriro/l) (see Eqs. (2.35), (2.39), and (3.48) and Figs. 3 and 10). These scaling 
functions are universal, i. e. independent of interaction potential parameters, and 
reflect the geometrical properties of the substrate (see Eqs. (A16), (B30), (C4), and 
(C5)). 

3. At liquid- vapor coexistence /i = Ho{T) and above the wetting temperature Tw \ro) 
(Sect. Ill A) the maximum fllm thickness Imax is determined by Eq. (3.2) depending 
on the dimensionless parameter k = {TQa'fg/aY^^ with the Hamaker constant a (see 
Fig. 4). In the hmiting cases k » 1 and k one obtains for l^ax power laws as 
given by Eqs. (3.9) - (3.11). 

4. For a small undersaturation A// = iiq{T)—ii > and above the wetting temperature 
Tw\tq) the equilibrium film thickness Iq on curved substrates equals the equilibrium 
film thickness on the corresponding planar substrate for the fiuid taken to be off 
coexistence at an effective undersaturation A/igjj. In the particular case of large 
radii rg, i. e., k — oo, the film thickness on a curved substrate equals that on the 
corresponding planar substrate as if the liquid-vapor coexistence curve was shifted 
by a constant value S^i'^ = ra^^ / (roAp) into the vapor phase (see Eqs. (3.9) and 
(3.28) and Fig. 5). In general, however, Iq and A/ieff are determined by k and the 
additional dimensionless parameter A = roApA/i/cr^^-' (Eq. (3.14)) which measures 
the actual distance A/i from liquid-vapor coexistence in terms of the difference 
Ap = pi — pg of the liquid and vapor number densities. This (k, A) parameter 
space separates into three regions (Fig. 6) in each of which A/^e// and thus the film 
thickness Iq reveals a characteristic power law behavior as function of k, and A (see 
Eqs. (3.29) - (3.31), (3.35), and (3.36) and Fig. 7). 

5. For suitably chosen interaction potential parameters the numerical analysis of the 
effective interface potential in its full form (Eqs. (2.10), (2.16), (2.17), and (2.29)) 
renders the phase diagram for first-order wetting of curved substrates (see Fig. 8(a)). 
The curvature of the substrate raises the wetting temperature T^\ro) compared to 
the wetting temperature Tif-* of the corresponding planar substrate. In the limit of 
large radii tq this shift vanishes as Tw\rQ) — Tw^ ~ ''"o At Tw\rQ) the jump A/q 
in the film thickness is finite (Fig. 8(b)) and diverges for ro — >^ oo as A/q ~ Tq . 

6. The curvature of a substrate erases a critical wetting transition such that there is 
a continuous and smooth increase of the film thickness (see Fig. 9). Even for very 
large radii tq the effects of curvature on the film thickness near this smeared out 
critical wetting transition are still significant (see Eq. (3.59) and Fig. 9(b)). At 
coexistence this film thickness is determined by the two dimensionless parameters 

/ , , N 1/3 

K = ( TQ(jig /a \ and 7 = 36/ (2aro) formed from the expansion coefficients a and h 
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of the planar effective interface potential (Eq. (2.27)). This (k, 7) parameter space 
separates into a rich structure (Fig. 11) such that its various regions correspond to 
different power laws for k (see Eqs. (3.59), (3.65), (3.68), (3.69), and (3.72) - (3.75)). 

7. Since the surface of a sphere and a cylinder is quasi-zero- and quasi-onedimensional, 
respectively, finite size effects smear out the first-order phase transitions shown in 
Fig. 8(b) as obtained by mean-field theory. However, according to the estimates in 
Sect. IV, these finite-size induced fluctuation effects are confined to a very narrow 
temperature interval 6T around the prcwctting line. For a typical spherical substrate 
with a radius Tq = lOOcr where a denotes the diameter of the fluid particles one flnds 
5T/T^^ ~ 10~^. For a cyhndrical substrate with the same radius 5T is vanishingly 
small. 



Acknowledgements 

We thank R. Evans, A. Hanke, and M. Napiorkowski for many helpful discussions. 



34 



APPENDIX A: EXPLICIT FORM OF THE EFFECTIVE INTERFACE 

POTENTIAL FOR A SPHERE 



According to Eqs. (2.17) and (2.29) a^^ and part of u;^^\l) are determined by t'^^^r; a) 
which represents the interaction potential of a fluid particle at distance r outside a sphere 
of radius a filled with the same fluid particles. In polar coordinates t^^\r; a) takes on the 
following form: 



&\r;a)= / d^r'wQr - r'\) 

J\r'\<a 



(Al) 



= 27rJ^ dr' r''^ dO sin 9 w Vr^ + r'^ - 2rr' cos 9^ , 



where we have already used the rotational invariance of t'^^\r] a). With the substitution 
§2 ^ ^2 _^ _ 2rr'cos9, Eq. (Al) yields 



a ^.1 t-r+r' 



&\r;a) ^ 2n I dr'- I ds sw{s) 

Jo ^ Jr-r' 

2^ rr+a POO 

— — / dr'{r — r') I ds sw{s). 

^ J 1 — a J r' 

Equation (2.3) leads to (on physical grounds only the case a > cr is of interest) 



(A2) 



tj)(r; a) = ^ JJ"^" dr'{r - r') dss(j)Lj{s), r>a + a 



~a So 
+ V Ir-adf'i'^ - f') la dsS(f)Lj{s), a < T < a + (7 



(A3) 



One has ty\a + cr, a) = t<''(a + a, a). In the next step we now explicitly use Eq. (2.4) so 
that 



a 



12 



1 



20 \r(r + a)^ r(r - a)^ 



2a 



12 



1 



+ 



45 \{r + af {r - af 

6/1 1 



a 



and 



2 \r(r + a)2 r(r — J 3 \(r + a)^ (r — a) 



t-^\r- a) = 7re -—(t + -ra^ - h -— + 



9 5 
1 1 
9 (r + a) 



5 r 2 r 



5 \8r(r + a)^ 



2eaV - 



1 1 



1 1 



2r(r + a)2 3 (r + a)^ 



(A4) 



(A5) 
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The integral J^°° dr t^'^\r; a), b > a, is given as 

Jl^^" dr rH\^\r; a) + (ir t>''(?"; a) , b < a + a. 



From Eq. ( |XB| ) one can directly infer the liquid-gas surface tension a^^^^ (Eq. (2.17)) 



(V^)2 (Vt^)H9 27 J 2073Q{h/a 



where a^^^ = ^ecr^(Ap)^. According to Eqs. (2.29) and (A6) one has to distinguish the 



'1,9 ~ 4 

cases I > don + a and d^ < I < d^ + a. For I > d^, + a one has 



io^'^ {l> d^ + a) 



(A^ 



whereas for < / < rf^ + cr 

cj^"*) {d^<l<d^ + a) 



Ap( pi r dr (-\ t^${r-r{)+pi [ dr ( -] t 



2 

>\r;ri) 



-Pv 



oo / ^ \ 2 



dr i^- ] v^'\r;ro) ]. (A9) 



If w'''^^(r;ro) is approximated as the linear superposition of (p^jir) (Eq. (2.5)) it has the 
same form as t>^(r;ro) with (e, a) replaced by {ewf,crwf)- For additional simplifications 
see Ref. [34(b)]. In this case and with Eq. (A4) we obtain 



{l>d^ + a) = Ap^piea^^ fl:'^l{ri, h) - pn^e^f(j]^ffl%{r(,, h) 

-piea'fiSM, h) + p^e^;a6^/it (ro, h)] (AlO) 

where 



and 



r(s\ / r N ^ fb'^ + Mb + a? b^ -8ab + a^ \ , . ^ , 

f''>-'>^ = i4^[-lbTaf (6315^1 (All) 



fl(a, '>) = :^( 2al,,^±4^ - In |±i^') , (A12) 
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Expansion in powers of — yields 



ro 



a;(^)(0=a;S(0+^il(0 (A13) 



where 



^ A 12 / 9^ ~ '^0 9/i + To ^ 



and 

In the hmit I ^ cr the functions Ss{ro/l) — ^attri^) I ^Mri^) — 
((1 + x-i) + f are given by 

Ss{x) = ^ (2x(x + If^^^^ - ln(2x + 1)) (A16) 

and 

1 1 + 2x 
2(x2(l + x))V3 ' 

respectively. 



APPENDIX B: THE EFFECTIVE INTERFACE POTENTIAL OF A 

CYLINDRICAL SUBSTRATE 

For a cylinder the calculation of the /-dependent terms a'fg and uj^'^\l) is somewhat 
more complicated than in the case of a sphere. They have the form 

<^S = -^(Ap)¥^n^,^;N) (Bi) 
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and 

u;^^\l) = Ap {pif^'\n, h- [w]) - p^g^''\ro, h; [v])) , (B2) 

with 



C(6,oo) C(0,a) 



and 



5^^Ha,&;M) = ^^ / / (B4) 

C(6,oo) C(0,a) 

where C(r<,r>) denotes the region between two cyhnders of radii r< and r>. Eq. (B3) 
can be rewritten in the form 

f^^^a, h- [w]) := I{a, oo) - /(a, b) (B5) 

where 

C(0,6) C(0,a) 

With r = (rj_, z) where the z-axis is the cyhnder axis and r_L orthogonal to it one has 
I{a, b) ^ dz J (Pr^ J (Pr'^w (^^{r± - r'^)' + z^^ 

0<|rx|<6 0<|r'j_|<a 

1 r°° 

=: / dzJ(z) (B7) 



with 

J{z)= J d^r^ J d\'^w{j{v^-r'^y + z^^ . (B8) 

0<|rx|<6 0<|r'_i_|<a 

The coordinate transformation 
yields 
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J{z) = 2nJ^ dffw (^Vf2 + ^2^ j ^2^_ (BIO) 

Equation (p9|) provides the following limits of integration for R: 

0<(^R+^f^ <6^ 0<(^R-^f^ <al (Bll) 

This corresponds to two circles, one with radius h centered at R = — the other with 
radius a centered at R = +|f. Thus f^d^R =: D{r,a,b) is the area of intersection 
of these two circles which depends on f for b — a < \r\ < b + a and equals vra^ for 
< |f| < 6-a: 

9a,b{^), b-a<\r\<b + a 
D{r,a,b) = { na"^, 0<|f|<6-a (B12) 
0, otherwise 

where ga,b{f^) := — ^^/—f'^ + 2n'^f'^ - + a? arcsin — b"^ arcsin "^^j^*"^ with rn? = 

52 — a2 and ^2 = 52 _|_ ^^2^ q^^^ j^^^g ^^^^ _ ^^2 i^gg^^gg ^j-^e smaller circle lies always within 

the larger one and 



with ga,a{f) = TTO^ — i^^/ ^0? — r"^ — 2o? arcsin ^. Inserting Eq. (P12|) into Eq. ( [B10|) yields 



J{z) = 27r^a^ / dssw{s) 

J\z\ 

V(6+a)2+22 



+2tx / dssw{s)ga,b [Vs^ -z^) (B14) 

J ^{b-ay^+z^ ' ^ ^ 



so that 



27ra 



"oo /" {b—a)^+z^ 

I{a, b) = / dz dssw{s) 

^0 Jo J\z\ 



and especially 



2 fV (fe+a)2+2;2 

+ — / (i^ / dssw{s)gafi iVs'^ — z'^) (B15) 



I{a, oo) = / dz dssw{s). (B16) 

'"o Jo J\z\ 
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Thus independent of the exphcit form of w f^^\a, b; [w]) is given by 



ro Jo 


/■oo 

/ ^ 










To Jo 





dssw{s)ga,b (y^^^) . (B17) 



Equation (2.3) requires the consideration of two cases: 

O / /»CXD /*00 

f^''\b> a + a) ^ -Ina' dz dssw{s) 

^0 V Jo J^/(b-a)2+z^ 



fOO p^/{b+^f+^ 



and 

/('=)(a <b<a + a) = 



/ dz dssw{s)ga,b(V^^) ) (B18) 

Jo JW<b-a)^+z^ ^ 'J 



27ra 



2 / j-^J<T'^-{b-a) 



ry <J^ — \o—a}' poo poo poo \ 

/ dz ds + dz / ds j sw{s) 

Jo Ju J Jc'^-ib-aY J Jib-aX^-^-/^ I 



^0 y/o J<j J ■yj c"^ -{b-aY Jy/{b-a)^+z 

— 1 / (i^; / ds + dz / ds | x 

'^0 Wo Ja J-sJa'^-ib-aY J {b-a^+z'^ 



X 



(y^^^) . (B19) 



For a = 6 one obtains 



rKa^b)^ 



2 

^0 VJo 



/ dz ds + dz ds \ sw{s) 

Jo Ja Ja Jz J 

poo p\ 

J a J z 

xsw(s)5„,„ (V^^^) . (B20) 



/•\/4a^+2^ poo p^/ToF+z^ 

dz I ds+ dz I X 

J a J a J z 



If we now adopt the Lennard- Jones form for w(r) (Eq. (2.4)) f'^^^ separates into a contri- 
bution frel from the repulsive part proportional to r~^^ and into a contribution f^ll^. due 
to the attractive part ~ r^^ of this potential so that 

f^'^ = ea'^m - eaVit (B21) 

and 
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with the coefficients gq = 8m^°, 02 = 72m^^n^, 04 = 36m^^ + 252m^^n^, oe = 252m^^n^ + 
420m^r^^a8 = GSm^^ + 630m^n^ + 315mV,aio = 70n^° + 620mV + 318m^n^ 012 = 
155nS + 55m^ + 462mV, cm = 132n^ + 156mV, aie = 52n'' + 20m^, Oig = Sn^ where 
= h"^ — a? and = 6^ + a^. By using a symbol- manipulation program one can show 
that the latter expression is identical with 



64ro69 " ^ V 2 2' ' V6 
The attractive contribution has a similar form: 

J attr - TTT -J \ . . . „ o o . :Ts79 ^ 



2ro Jo V (z^ + 271^^2 + m^f^ 

2ro632^^ 1^2' 2' ' b\ 

For a = ri and 6 = /i one obtains f^'^^ri, h; [w])- 

In order to provide an explicit formula for the effective interface potential lu^'^\1) one fi- 
nally has to specify the contribution g'^'^^ from the substrate potential (sec Eqs. (B2)) 
and (B4). If the substrate potential is approximated as the linear superposition of 
Lennard- Jones interaction potentials between pairs of substrate atoms and fluid atoms 
(see Eq. (2.5)), g^'^^ has the same functional as f^'^^ with e and a replaced by e^/ and a^f, 
respectively: 

g^^'^r-o, h) = /(^) (a = ro, 6 = /i; e ^ e,/, a ^ <j^f) . (B25) 

Equations (B23) and (B24) are only valid for h > tq + rf^ + a, because due to the 
Heaviside function entering Eq. (2.3) f'^^^ changes its form for d^ < I < d^ + (J- Since 
in the substrate potential there is no such Heaviside function, g^'^^ retains its form for all 
values of h. Keeping this in mind one has 



>d^ + a)^Ap {pif^^^ (n, h) - p^f^'^^ (ro, h; e ^ e^f, a ^ a^f)) . (B26) 



In the limit ^ — >■ Eq. (B26) yields exphcitly 



= + a;ll(/), d^/ro « 1, (B27) 



where 
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^^'A...^i2ro p fll 9... /ro\2\ 99r2 /13 11 /ro\2\\ rf. 



I ( ^ ) I (B28) 



and 



^attr 



^\ 6^0 p /5 3 /ro\2\ 15 ^ /7 5 /ro\2\\ 
-T^^^'^"/^i2 2F,^-,-;2; (-) j + T/I^ ^F, -; 3; (- j jj- 

The functions 5c(3:^) and Sc{x) (Eqs. (2.36) and (3.4)) are thus given by 

3n X „/53„ / 



and 



1/3 



/4(i+x)^y 



respectively. 

The hquid- vapor surface tension follows from Eqs. (Bl) and (B20). A lengthy calcula- 
tion yields the following result based on Eqs. (2.3) and (2.4): 

= -{^pfea'K {^^ (B32) 

where 

4 

K{x) = ^Ki{x) (B33) 

with 
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Ki{x) 



+ 



TT V4x2 + 1 



35 



125 



337 



209 



+ 



2048a; 49152a;3 786432x5 2097152a;7 
1 3 \ 



+ 



— TT 



1024x3(4x2 + 1)3 12288x5(4x2 + 1)2 16384x^(4x2 + 1). 
175 In (2x + V4x2 + l) 



4194304 



x"- 



(B34) 



K2{x) — J dz<^ 



( 35x6 ^ 3Q^4^2 ^ g^2^4 ^ ^6 



^9(4x2 + ^2)7/2 



X 



X I I - arcsin 1.1(1- z^) (l + ^] ) ) + arcsin 



+ ^V(l--^)(4x2 + .2-l)(l - _L 



14x2 + 3^2 



12z'i(4x2 + ^2) 
70x^ + 32x2.^2 ^ 420a;6 ^ 290x^2 + 62x2^2^ + 3^^ 



12^6(4x2 + ^2)2 



12^8(4x2 + z'^Y 



)}■ 



(B35) 



and 



( \ ^ rr^i^ ^ 1 V4x2 + 1 + 2x 

-f^sl^j = V 4x"' + 1 — -— ^ m , , 

^ ^ 8x 16x2 V4x2 + l-2x 



(B36) 



K^{x) — 2x J dz arcsin 



1 ITT 

+ 



^r^\ 1 (1 - .22)3/2^4x2 + Z^-1 



2x / 4x2 z^ 

{2x^ + z^){l-z^)-2x^z^ 



— arcsin 



2^3^4x2 + ^2 \^2 2x2 
In the asymptotic limit /i — > oo this expression reduces to 

<7g(/.»a) = 
(,)/ lln(V^_^/5 ^ 
V 6 {h/ay [h/af VS 144 

^ 1 +0^ ^ 



1536 (/i/t^)' 



{h/a) 



(B37) 



(B38) 
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APPENDIX C: EFFECTIVE INTERFACE POTENTIAL FOR COATED 

SPHERES AND CYLINDERS 



We consider spherical and cylindrical substrates which are composed of two different 
species for < r < tq— ^ and tq—S < r < tq. They interact with the fluid particles accord- 
ing to pair potentials and respectively. For reasons of simplicity we assume 
that the number density of the substrate particles is constant. If one approximates 
the substrate potential y(r;ro) = Pu,i^(r;ro) by the homogeneous (i. e. disregarding the 
discrete arrangement of the substrate atoms) linear superposition of these pair potentials 
one has v(r; tq) — v{r] Tq) with 



v{r;ro) = f dV0i°j(|r - r'|) + / 

^5(r-o-<5) ^5(ro)\5(ro-<5) 

= / ^V02(|r - r'l) + / U:j{\r - r'|) - 0«(|r - r'|)) 



5(ro) JS{ra-5) 

=: T;i(r;ro) +T;o,i(r;ro - 5) (CI) 

where iS(r) denotes a substrate of radius r. If the functions 0^]-(r) and 0|^]-(r) are ap- 
proximated by Eq. (2.5) with interaction potential parameters and o"!)/}, 
respectively, the computation of Wi(r;ro) and fo,i(r;ro — S) can be carried out along the 
hues of the Appendices A and B. On this basis we obtain 

-P. {^t]{<yt]f - ^S(^S)') &Xr, -S,h)+ p.6S(a«) Vi?*?(ro, h) 
+piea''fj:^,f{n, h) - piea'fitC\n, h)) . (C2) 

The expansions of the effective interface potentials cj^'^^l,ro) into powers of Tq/I under 
the condition that both 5 <^ ro,l and du, <^ ro,l can be calculated straightforwardly. As 
a generalization of Eq. (2.27) one obtains 

where Sr{x) and TZr{x) are given by Eqs. (2.36) and (3.51), respectively, and with the 
constants a and h given by Eqs. (3.44), (3.47), and (3.48); they are identical to those 
of the corresponding planar substrate. For the model potentials given by Eqs. (2.4) and 
(2.5) one has 

^/vStt 1 ^/53^/x,, 



2 + ' I 2'2' ' \x + l 
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for a cylinder and 



8{l + xf 
(1 + 2x)3 



(C5) 



for a sphere, respectively. From these results one can infer the explicit form of the functions 
qr{x,'y) (Eqs. (3.55) and (3.56)): 



<lc{x) 



3x 



X 



^.53 
2F1 -,-,2, 



2' 2' ' \x + l 



X 



5 a;2 I. / 7 5 



4(x + l)2 



2' 2' ' \x + l 



X 



X 



^.53^ 
2F1 I ^,^,2 



2' 2' ' 1 



X 



8 (x + l)2 



2' 2' ' + 1 



a; 



+ 



175 i- i I ^ 7 f X 

+^^— TT2F1I -,-,4 



48 (x + 1)' 



(C6) 



for a cylinder and 



1 1 + 2a; + 2x2 
^^("^ = x 1 + 2X 



(C7) 



for a sphere, respectively. 
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Figure captions 



FIG. 1. (a) Cross section of the density configuration around a sphere or a cyhnder of 
radius tq with a number density p^. The gas phase with number density Pg is truncated 
for r > L. The hquidhke film with number density pi has a radius of curvature given by h 
and a thickness I = h — tq. The density distribution vanishes within an excluded volume 
To < r < ri = rQ + dyj. (b) Corresponding density configuration close to a planar wall. 

FIG. 2. Liquid- vapor surface tension of a liquid drop (s) and thread (c) as compared with 
the planar surface tension (p) as function of their radius h in units of the diameter a of the 
fluid particles based on Eqs. (2.16) and (2.17) for the model potential given in Eq. (2.4) 
with the parameters cr = 1.0 and e = 1.0. Within the sharp-kink approximation the ratios 
o'l^g (h) / al^g are constant as function of the temperature and behave asymptotically for 

a/h as given in Eqs. (2.22) and (2.23). For a/h < 1.89-10-^ one has (^I'^^ih) < a{J(/i); 

otherwise <j\^g{h) is smaller than crl^g{h). Beyond the sharp-kink approximation the ratio 

<^\^g{h)/a'fg for a sphere is expected to approach 1 linearly as function of cr//i — ^ with a 
slope —25/(7 given by the Tolman length 5. (See also Ref. ||59[| .) 

FIG. 3. (a) Scaling function St-{x) = iS^j2-(^) of the effective interface potential (see 
Eqs. (2.29), (2.30), (2.32), (2.34), (2.36), and (2.39)) as function of x = tq//. S^ix) 
vanishes linearly for x ^ (Eq. (2.41)) and approaches 1 for x ^ oo proportional to 
(Eq. (2.40)). The maxima (•) and the inflection points (x) are located at Xmax — 
2.56, Xmax — 1.85, and x^f* ^ 4.24, x^^'' ^ 3.18, respectively, (b) Corresponding scaling 
functions Sr{x) (see Eq. (3.41)). According to this figure Eq. (3.3) has indeed one unique 
solution. 

FIG. 4. Ratios Imax/f'o = V^* ^ ^ function of /t = {rQcrl^g /a)^/^ obtained by solving 
Eq. (3.3) numerically for the same interaction potential parameters as those used in Fig. 2. 
According to Eqs. (3.6) and (3.7) both ratios vanish for k — > oo proportional to 
whereas for /t — > they diverge proportional to k"^/^ and for the cylinder and the 

sphere, respectively. 

FIG. 5. Complete wetting film thickness along an isotherm ksT/e = 0.94 for a sphere and 
cylinder, respectively, of radius lOOOo", compared with that of a planar wall (ro = oo). 
The curves correspond to Eq. (3.14) with the interaction parameters chosen to be the 
same as in Fig. 2. In this case n ~ 75.5 so that Eqs. (3.21), (3.25), and (3.26) apply. 
Therefore the complete wetting film thickness on a curved substrate resembles that on the 
corresponding planar substrate as if liquid-vapor coexistence would occur at A/i(K, A) = 
—6p{K,X = — r) < (see Eqs. (3.17) and (3.27)) instead of at A/i = 0. (A/ig// = 
implies A = -r, see Eqs. (3.24)-(3.26).) 

FIG. 6. According to Eqs. (3.25) and (3.26) the thickness of a complete wetting film on 
a curved substrate is given by the film thickness on the corresponding planar substrate 
for an effective undersaturation A/ie//(A/U; tq, r). In terms of the relevant dimensionless 
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parameters k = {(y^^gT^/ a)^/^ and A = AproAfi/al^g one finds three distinct regimes I, II, 
and III within which \eff{n, A) exhibits the characteristic power laws given by Eqs. (3.29), 
(3.30), (3.31), (3.35), and (3.36), respectively, (a) shows the crossover lines kjji, khjh, 
and Kjjii between these regimes for a cylinder, (b) the corresponding ones for a sphere. 
Their amplitudes are chosen such that they meet in a single point. The variation of 
Ae//(/t, A) along the dashed-dottcd paths is shown in Fig. 7 

FIG. 7. (a) shows the effective undersaturation Ae//(/t, A) for a cylinder along the dashcd- 
dotted path in Fig. 6(a), i. e., for A = 0.35, and (b) for a sphere along the dashed-dotted 
path in Fig. 6(b), i. e., for A = 0.5, as function of k,. The full curves correspond to the 
full solution determined by Eq. (3.26). These solutions display the crossover behaviors 
which occur according to the corresponding vertical paths in Fig. 6. For large values of k 
one is in regime II where A^jj ~ r + A is constant as function of k (dashed curves). Upon 
decreasing k. one enters regime III where Eqs. (3.35) and (3.36) (dashed-dotted curves) 
hold. Ultimately for k ^ one reaches regime I where Eqs. (3.29) and (3.30) (dotted 
curves) are valid. Ag// is raised to the powers —5/6 and —2/3, respectively, so that the 
dotted curves are straight lines. 

FIG. 8. (a) Phase diagram of a first-order wetting transition on a cyhndrical and a 

spherical substrate of radii ro/a = 100 obtained by a numerical minimization of Eq. (2.10) 
with the substrate potential parameters cr^f/a = 1.35, e^f/e = 0.625, d^/a = 1.175, and 
PvjC^ = 1.0. The prewetting line of the corresponding planar substrate (solid curve) joins 
the liquid- vapor coexistence curve /j, — Ho{T) = tangentially at the first-order transition 
temperature T^^. The wetting behavior on a curved substrate is similar to that of a 
planar substrate whose coexistence curve is shifted upwards. Consequently for a curved 
substrate at coexistence there is a thin-thick transition at Tw\k) where the remaining 
part of the prewetting line intersects the coexistence curve linearly. Note, however, that 
along the prewetting line of the curved substrates k varies as function of the temperature 
and is thus not a constant; for the given parameter values k ranges between 13.8 and 17.1 
for 0.78 < hsT je < 1.0. Therefore the prewetting lines belonging to the curved surfaces 
are not parallel to the corresponding prewetting line of the fiat substrate. The thin- 
thick transition can be traced into the oversaturated vapor region (dashed-dotted lines) 
provided it is possible to maintain the thermodynamically instable vapor phase for these 
values of At the upper end points of these curves the minimum in Jls(/) corresponding 
to the thicker film ceases to exist. Each prewetting line ends at its lower end in a critical 
point. However, in the present model and approximation the prewetting lines of both 
the planar and the curved substrates can only be followed up to the dotted curve which 
denotes the loci where in the bulk the metastable liquid phase ceased to exist. Therefore 
beyond the dotted curve the present sharp-kink approximation (Eq. (2.6)) for the density 
distribution close to the wall is no longer applicable. In that region of the phase diagram 
for the present choice of interaction potentials the thicknesses of the wetting films are so 
small that these films cannot be described adequately by a liquidlikc layer with a number 
density pi which corresponds to that of the metastable bulk liquid phase, (b) Temperature 
dependence of the equilibrium film thickness Iq at coexistence for the systems described 
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by the phase diagram in (a). It illustrates the shift of the wetting transition temperature 
and the finite discontinuity induced by the curvature of the substrate. 

FIG. 9. (a) Schematic phase diagram for a planar substrate exhibiting a second-order 
wetting transition at Tw\ At /i = i^o{T) and for T < Tc the bulk hquid and gas phases 
coexist; only temperatures above the triple point are considered. If one records the 
film thickness Z/a as function of the temperature ksT/e along a path slightly below 
coexistence (thin line) and at a fixed undersaturation > (dotted line) for a given 
set of interaction parameters (see below) one obtains the curves shown in (b). For the 
film thickness on a cylinder or a sphere the path along the coexistence curve across Tw^ 
is equivalent to the dotted path for the corresponding planar substrate because of the 
effective upward shift of the bulk phase diagram due to the curvature (see main text), 
(b) Film thickness //o" as function of the temperature ksT/e for a planar substrate and 
for a cylinder and a sphere of radius r^/a = 10^, respectively, on a thermodynamic path 
along the coexistence line, = fio(T), crossing the second order wetting transition of the 
planar substrate. The interaction parameters have been chosen such that for the planar 
substrate there is a wetting transition at ksTw^ /e = 0.9. This requires (Eqs. (3.41) - 

(3.46)) ei°}/e = (a/ai°]) V(rirVp«. - 0-6002 and a'^j/a = 1.0 so that a (t = xi"^^ = 
(Eq. (3.41)) is fulfilled. In order to satisfy the remaining conditions for a critical wetting 
transition (Eq. (3.41)) one may choose, e. g., e^|/e = 1-4, o'^j/'^ = 1-0, and 6/a = 1.5. 
The inflection points (•) in the film thickness l^'^\T) on the curved substrates occur below 
the wetting temperature Tw^ of the planar substrate. The levelling off of the curves for 
the cylinder and the sphere is not yet complete even at Tc. It is remarkable that for 
)/Tif) ^ 10-2 even for radii of the order of 1cm the effect of the curvature 

on the wetting film thickness as compared with a truly planar geometry is still visible. 
We have checked that the above choice of interaction potential parameters is not atypical 
in the sense that they cause an unusually small temperature gradient of the Hamaker 
constant at Ty^ as compared with the experimental observation of critical wetting |65 |. 
Therefore we conclude that this extraordinarily strong sensitivity of critical wetting to 
curvature is a generic feature. 

FIG. 10. Scaling function 7lr{x) given by Eqs. (3.49), (C4), and (C5), as function 
of X = ro/l. TZr{x) rises linearly to a finite value for x — {lZc{x = 0) = ^ and 
TZs{x = 0) = 8, respectively; Eqs. (3.51) and (3.52)) and approaches 1 for x oo 
proportional to x^^ (Eq. (3.50)). 

FIG. 11 In the case of critical wetting at coexistence the film thickness on the curved sub- 
strates is determined by Eq. (3.54) with the dimensionless parameters k = [rQCTig /a) 
and 7 = 36/ (2aro). Since x = tq/I is positive, it is sufficient to consider the regions (k > 0, 
7 > 0) and (k < 0, 7 < 0) which correspond to T > T^f'' and T < Tw \ respectively. The 
limits T = Tw^ ± correspond to 7 — ±00 whereas on the line 7 = one has, apart from 
possible multicritical wetting transitions, first-order wetting which we have studied in 
Subsect. inc. By means of the limits x ^ and x — >■ cx) of the functions qr{x) and s^(x) 
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entering into Eq. (3.55) one can evaluate Eq. (3.54) for the regions I, II, and III of the 
(k, 7) parameter space. In the regions Ia,b one has :§> 1 and I7I » 1 so that 7 <S k^, 

i. e., T — Tif"* ± and Tq ^ ^36/ (2o"[^^) j , which corresponds to film thicknesses given 

by Eqs. (3.57) - (3.59). For \k\ ^ 1 and I7I ^ 1 with k'j <^ 1 (large Tq and 6/a — 0) 
the film thicknesses are determined by Eqs. (3.60)-(3.62) for positive 7 (region Ila) and 
by Eqs. (3.64)-(3.65) for negative 7 (region lib). In the case \k,\ <S 1 (regions IIIa,b), i. e. 

ro <^ (^\a\/aj^g^ , Eqs. (3.66) - (3.69) are valid. For I7I — > 00 (regions IIIai,bi), i. e. 

T — > Tif^ ± 0, one can expand Eqs. (3.66) and (3.67) into powers of K^/7 2rQal^g /{3b) 
which leads to Eqs. (3.70) - (3.73). The opposite hmit, 7 — > makes sense only for 7 > 
(region IIIa2) where the film thicknesses are determined by Eqs. (3.74) and (3.75). In the 
remaining regions IVa,b and Va,b there is no suitable expansion parameter so that there 
Eq. (3.54) can only be solved numerically. 
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